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^ ! Abstract 

The eigenvalue problem of the spherically symmetric oscillator Hamiltonian is revis- 
' ited in the context of canonical raising and lowering operators. The Hamiltonian is 

^-j. ■ then factorized in terms of two not mutually adjoint factorizing operators which, in 

CN . turn, give rise to a non-Hermitian radial Hamiltonian. The set of eigenvalues of this 

new Hamiltonian is exactly the same as the energy spectrum of the radial oscillator 
■ and the new square-integrable eigenfunctions are complex Darboux-deformations of 

the associated Laguerre polynomials. 



q-i 1 Introduction 

Factorization is a powerful method to compute eigenvalues and recurrence relations for 
CO ' solutions of differential equations in mathematical physics (see e.g. [1-5]). Introduced 

in Quantum Mechanics (QM) since Fock, Schrodinger and Dirac times, some of the an- 
tecedents of this algebraic method can be found in the geometrical formalism of Darboux 
\ and Backhand. It is also notable that the (canonical) Fock-Dirac ladder operators have 

been of big influence in the development of contemporary Physics over the years [6] . The 
method, as formalized by Infeld and Hull [7] , was extended by Mielnik to embrace factor- 
izing operators of the harmonic oscillator Hamiltoninan which do not create nor annihilate 
occupation numbers in the Fock states [8]. Applications of the Mielnik's method to radial 
problems were reported first by Fernandez [9] and connections with the Supersymmetric 
formalism of QM were found by Nieto [10] and Sukumar [11]. The relationship between 
factorizing and Darboux-deforming potentials in QM was extensively studied by Andri- 
anov and coworkers [12]. Nowadays, the factorization method is of the highest utility 
in constructing new exactly solvable potentials in QM (see e.g. [13, 14] and references 
quoted therein). In particular, recent results show that the method can be used in solving 
the problem of complex potentials with real spectrum [15-20] as well as in the study of 
resonances and Gamow-Siegert states [21-23]. 

In the one-dimensional case there is a unique expression for the canonical ladder oper- 
ators since the energy levels are not degenerated. For instance, the Fock-Dirac operators 
transform an eigenfunction of the harmonic oscillator into another one and factorize, at 
the same time, the corresponding Hamiltonian. However, in a general situation, the pair 
of operators which factorize a given Hamiltonian do not necessarily correspond to the 
shift operators. Distorted versions of the harmonic oscillator, for example, require higher 
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order raising and lowering operators although the Hamiltonian is factorized by first order 
differential operators [8,24]. On the other hand, it is usual to find degenerate energy levels 
in higher dimensions and several kinds of shift operators are necessary to intertwine the 
wave-functions [2]. A particular case is of special interest: spherical symmetry in central 
potentials induces energy degeneracies under rotations (different orientations of the angu- 
lar momentum vector lead to the same energy). Hence, the ladder operators are labelled 
by the azimuthal quantum number £ and their action mainly affects the ^-dependence of 
the solutions (see for instance [7,9,25-27]). Thus, they intertwine solutions of potentials 
which differ in one unit of the angular momentum. 

In this work we are interested in the three dimensional isotropic oscillator. The Hamil- 
tonian of this system can be factorized in four different forms by means of a basic set of 
canonical ladder operators [28]. Each of these products preserves the adjointness of the 
Hamiltonian since the factors are mutually adjoint. In contrast, we are going to express 
the radial oscillator Hamiltonian as the product of two not mutually adjoint factorizing 
operators. The advantage of these new operators is that they give rise to non-Hermitian 
radial Hamiltonians for which the point spectrum is exactly the same as the spherically 
symmetric oscillator one. 

In Section 2 conventional factorizations of the isotropic oscillator are revisited. Ex- 
plicitly derived, the four canonical ladder operators are used to construct the basis of 
physical wave-functions (orthonormal, associated Laguerre polynomials) and the energy 
spectrum. For a given £, the creation or annihilation of nodes in the wave-functions 
require second order differential operators which do not factorize the Hamiltonian. In 
Section 3 a pair of not mutually adjoint factorizing operators is introduced to extend 
the number of factorizations in which the radial oscillator Hamiltonian can be expressed. 
These same operators are used to construct a new Hamiltonian which is not self adjoint 
under the inner product of the Laguerre-Hilbert Space H. The relevant aspect of our 
results is that the eigenfunctions of the new Hamiltonian are square-integrable in Ti. It 
is also introduced a pair of ladder operators which increase or decrease the energy of this 
new system. As we shall see, these shift operators are necessarily fourth order differential 
operators. Finally, the paper is closed with some concluding remarks in Section 4. 

2 Radial Oscillator Revisited 

The time-independent Schrodinger equation for the isotropic oscillator V(r) = r 2 , in 
appropriate units and after separation of variables, reduces to 



where the azimuthal quantum number £ is a non- negative integer, E(£) is twice the di- 
mensionless energy eigenvalue and the radial Hamiltonian reads: 




(1) 



H e = 



dr 2 



+ V e (r) 



dr 2 



+ 



1(1+1) 



(2) 
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The effective potential Vi(r) has the domain Dy = [0, oo) and the functions <j>(r,£) are 
defined in terms of the usual radial wave-function <f>(r,£) = rR(r,£). To be physically 
interpretable, these last functions have to satisfy J +o ° \R(r, £)\ 2 r 2 dr < oo. In the sequel, 
whenever there is no confusion we shall use the shortcut notation / = f(r, £, Ef, . . .), keep- 
ing implicit the dependence of / on r, £, ee, and other possible variables and parameters. 

Conventional factorization is useful in solving equation p]) by expressing (j2J) as the 
product of two mutually adjoint first order differential operators 

H e = ala £ + e e ; a\ : = + a(r, £); (a]) 1 " = a t (3) 

with Ee a real (factorization) constant to be fixed and a(r,£) a function satisfying the 
Riccati equation 

- a ' + a 2 = V e - e e . (4) 
A particular solution of equation is immediate: 

a = r-*-^, e e = 2£ + 3. (5) 
r 

This superpotential a leads to one of the canonical forms of the factorizing operators: 

d £+1 t d £+1 

at = h r, a\ = — ; h r. (6) 

dr r dr r 

The first advantage of the method is clear by noticing that a solution of a£(j)o{f,£) = 
gives a solution of <^f belonging to the eigenvalue E (£) := E$ = Eg. The calculation 
leads to 

MrJ):=4\r) = Ci ) £} r^e- r2 / 2 (7) 

which is free of nodes and becomes zero at the edges of Dy. Hence, Cg can be taken 
as the normalization constant and, by applying the Sturm oscillation theorem [29], we 
know that Eq = 2£ + 3 is the ground state energy of H^. Now, it is useful to reverse the 
product to get 

H l+ i = a e a\ + ee - 2 := a<> a\ + $ e . (8) 

Notice that the azimuthal quantum number £ plays the role of a parameter in all the 
previous expressions. Hence, by changing £ — > £ — 1 and £ — > £ + 1 in (jS]) and ([3]) 
respectively we also obtain 

— a,e—\0'i_i + Hp+i — a e+i a v+i "I" ££+!■ (9) 

As a consequence, the following intertwining relationships are true 

a t (H t -2) = H t+1 a t , a{_ x {H e + 2) = H t . x a\_ v (10) 

Thereby, if <f>(r, £) is eigenfunction of Hi with eigenvalue E(£), then <f>(r, £ + 1) oc 0(r, £) 
and <p(r, £ — l)(x a\_ x <p(r, £) respectively satisfy the eigenvalue equation of i^+i and Hi-\ 
with eigenvalue E(£ + 1) = E(£) - 2 and E(i - 1) = E(£) + 2. 
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Let us take full advantage of these last results. After applying the R.H.S. relationship 
flTUJ) on the ground state <j>\ , , we arrive at the function 

<P ( t l) oc a\_ x 4>f = (-2r 2 + 2£ + 1) (11) 



which has a single node at r = ^/ [2£ + l)/2 e Dy and satisfies the Schrodinger equation 

H i - 1 *t 1) = (E!P+2)tf- 1 \ (12) 
Function (ITT]) can be rewritten as 

<t>t 1] oc ^ jr- 2 ^ 1 ^) ^ f _LjT) r 2(^i/2) e - 2 J . (13) 

After an intermediary change of independent variable r 2 = x, one arrives at 

0f- 1} (r) = Cf- 1} <pt l \r) Lf~ 1)+1/2 (r 2 ) (14) 

where C*j £ ^ is a normalization constant and 

1 r/ n 

jZ( x ) = — x- v e x —(e- x x n+u ), v > -1 (15) 
is the associated Laguerre polynomial of degree n. A similar procedure leads to 

4- 2 \r) = Ct 2) 4~ 2 \r) Lt 2)+1/2 (r>) oc a\_ 2 ^\r) (16) 

which is the eigenfunction of Hi~i with eigenvalue E 2 ^ = e[ 1 + 2 = Eq + 4. The 
procedure can be applied s times to get the eigenfunction of H e _ s with eigenvalue Eg S ' = 
E^ + 2s: 



's-1 



</T s) = Ct s) I! # : = C ^ S) { a ls ■ ■ ■ 4-2 <4-i) \ * = 1, 2, • ■ • (17) 

\fe=o / 

Changing £ £ + s and using Lq = 1, this last expression leads to 

tfM^^WW), , = o,i,... (is) 

Function (ITS]) is the physical solution of (jTJ) with eigenvalue = Eq +s ^+2s = E£ +s +2s = 
2(2s + £) + 3. The presence of nodes in <f>y is due to the zeros of the polynomial Ll +1 ^ 2 } 
the number of which is equal to s. By taking n = 2s + £ = 1,0, . . ., we can use the 
conventional notation <f>P(r) = 4> n (f) and Eg = E n — 2n + 3. 

Remark that the action of ai on an arbitrary function (jl8]) . different from (j) , produces 
£i£ = <Ps-\* ■ Thus annihilates one node of the wave-function at the cost of increasing 
in one unit the azimuthal quantum number £. This can be represented in the Cartesian 
(s x £)-plane as the diagonal displacement (s,£) — > (s — 1,£ + 1); a schematic picture 
of which is shown in Figure UJ The action of a\__ 1 reverses the previous result; this 
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Figure 1: Lattice representation in the (s x £)-plane for the action of the four canonical 
factorization operators ag, a\, be, b\, of the radial oscillator. 



is represented by (s,£) — > (s + 1,£ — 1). Notice that a particular situation occurs for 
£ = 0: The R.H.S. relationship (1101) prohibits the action of aj on any eigenfunction (f)^ 
of H because it produces the forbidden value £ = — 1. However, function ffT8]) can be 
also derived by applying the appropriate product of annihilation operators {ag} on a 
convenient excited state of Hq. The straightforward calculation leads to 

#=C?(no/-i-*) A ^0, 5 = 0,1,... (19) 
Wo / 

Another canonical factorization, different from ([3]) and (Q, is still possible for Hf. 

H e = hbl + K £ , bl:=~+j(r,e), (b e y = b e (20) 
where the factorization constant Kg is to be fixed and the superpotential 7(7", £) fulfills 

i + l 2 = V t - K t . (21) 
The immediate solution of this last equation reads 

7 = r + — , K e = -e e = -(2£ + 3). (22) 
r 

Hence, the reverse product (1201) is given by 

H e+1 = b\b e - Bi (23) 
and the new intertwining relationships are 

b^ (H e - 2) = H t -x 6/-1, 6] (#< + 2) = H i+1 b\. (24) 
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The action of b\ on the eigenfunction cpl of He with eigenvalue E s produces a square- 
integrable eigen-solution of He+i belonging to the eigenvalue Eg + 2 = 2(2s + 1 + 1) + 3 = 
Es . Notice that the number of nodes is preserved; we write (pi = Ci £ b\ <$p . Such 
an operation is represented in Figured] as the horizontal displacement (s,£) — > (s, £ + 1). 
On the other hand, the function <pi~^ = Cf ^ 0s is eigenfunction of He-i with 
eigenvalue Ef l ' = E^ — 2. Thus, the action of be~\ is represented by the displacement 
(s, £) — »• (s, £ — 1) in the (s x £)-plane, as it is also shown in Figure [Q 

Expressions (J3J) , ([9]) , (120]) and (123]) recover the canonical factorizations of reported 
by Fernandez, Negro and del Olmo in 1996 [28] (see also [30, 31] and references quoted 
therein). As we have shown, the four canonical ladder operators at, a], be, b\, induce 
displacements in the (s x £)-plane by either creating (annihilating) a node in the wave- 
function at the cost of decreasing (increasing) the azimuthal quantum number in one unit 
or by preserving the number of nodes but changing the value of £ in one unit. Although 
all these transformations involve a horizontal component in the displacements of the point 
(s,£), purely vertical displacements are also possible. Indeed, we have a pair of mutually 
adjoint second order differential operators defined as follows: 

(Sty = S e := ae-i b^ x = be a t 

(25) 

= -He + 2r£ + 2r 2 + 1. 

These operators increase or decrease in one unit the number of nodes in the wave-functions 

Sj^octfi, fl^Woc^ (26) 

and intertwine eigenfunctions of He for which the energy differs in four units 

S t He = (He + A)S e , Sj H e = (H e - A)Sj. (27) 

Hence, we get the following commutation rules 

[S t , He] = ASe, [S\, He] = -AS\. (28) 

It is clear that the Hamiltonian He is not factorized by St and S\: 

Se S\ = (He -s e + A) (He + e t ) , Sj S t = (H e +e t -A) (H e - e t ) . (29) 

However, there is a simple commutation rule 

[St, S}] = 8H e . (30) 

Figure [2] shows the schematic representation of two arbitrary but contiguous energy levels 
of He; the corresponding points of the (s x £)-plane are also depicted. As indicated, the 
operator Sj scales up the energy to the closer excited one and Sj operates in the opposite 
direction. 
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Figure 2: The ladder operators S\ and Si respectively create and annihilate nodes in the 
eigenfunctions of He. Notice that displacements in the (s x £)-plane are purely vertical, 
thus the value of the azimuthal quantum number £ is fixed. 



3 Complex Factorization 



In this section we are going to go a step further in the factorization of Hamiltonian 
Let us analyze the product 

H i = A i B t + e t (31) 

where ee is now a complex factorization constant with non-trivial imaginary part (from 
now on 9?(z) and ^s{z) respectively stand for the real and imaginary parts of z G C while 
z corresponds to its complex conjugate). The new factorizing operators are established 
in the context of a 'refined factorization' [32] to read 

At = ~+f3(r,i), B e = ^- + (3(r,£) (32) 
ar ar 

with (3 a complex function satisfying the Riccati equation 

- + (3 2 = V e - ee. (33) 

The first aspect which distinguishes (|3ip from is that Ae and Be are not mutually 
adjoint in the Hilbert space TC = L 2 (Dy) spanned by the vectors (fT4l) . However, as He is 
selfadjoint, He = B\ A\ + ee must be true. The straightforward calculation shows that this 
last expression leads to the complex conjugation of equation ( 1331) . the solution of which 
is immediate since Ve(r) is real. 

As in the previous section, the factors in (1311) lead to solutions of the involved eigen- 
value problem. In particular, the function u(r,£), annihilated by Be, is also solution of 
the Schrodinger equation 

Hiu = ee u. (34) 

Since ee is complex and He is an Hermitian operator, it is clear that u is not square- 
integrable in Dy. However, this function is very useful because the non- linear Ric- 
cati equation (1331) can be mapped into ( 1341) by means of the logarithmic derivative 



(3 = — -r-lnw. The appropriate transformations allow to express the solutions of (J3 
in terms of confluent hypergeometric functions i-Fi(a, c;z) (see e.g. [3,4]). In particular 
we shall use 



u(r, i) = c e r^e-^l \F, [i + \ - p + | r 2 
with C£ an arbitrary integration constant which will be fixed as 1. Hence, we have 



P( r ,£) = r- 



1 



2r 



2£ + 3-e e \ xF x (f + \ - f ,1 + f; r 2 ) 



4^ + 6 ; ^(f + S_a i+| ;r 2) 
The well known analytical properties of the \Fi functions lead to 



(35) 



(36) 







(37) 



-oo 



The reverse product fl3T|) produces a new second order differential operator 

B e A e + e e = H e + 2f3> (38) 
which is a non-Hermitian Hamiltonian in TC since /3 is complex. The new potential 

v(r,£,e £ ) :=V e (r)+2f3'(r,£) (39) 
is mainly a real function at the edges of Dy\ 

v (r^0,£,e e )=V i+1 (r^0)+i^ 

v(r —> oo, £, e<?) = \4(r ^ oo) - 2 



(40) 
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Figure 3: Real (left) and imaginary (right) parts of the complex potential fl39|) obtained 
from Vo(r) with a factorization constant e = 3 + i 10~ 3 (blue curve), e = 7 + i | (red 
curve) and e = 11 + i5 (black curve). Observe that greater values of 9ft(e ) induce larger 
displacements to the left in the deformation 



Figure [3] shows the global aspect of the new potentials v(r, £, eg) for I = and several 
values of 6q. Notice that the Darboux-deformation (1391) of the initial potential Vg is 
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characterized by a zone M C Dy surrounding the origin. Namely, M = [0,tq), with tq a 
finite number depending on i and e^. In this deformation zone M, the function dt(v) is 
merely a distortion of V^+i while 3(f) is a nontrivial position-dependent function. For the 
sake of notation, we shall write V£ + x for the complex potential defined in (I39p . Thereby, 
the non-Hermitian Hamiltonian (|38p reads 

d 2 

h i+1 = B t A t + e e = --7-5 + v i+1 . (41) 

Once we have obtained the factorizations (13"Tj) and (|4ip . it is immediate to arrive at the 
following intertwining relationships: 

Be Hi = h e+1 B e , An h e+1 = H e A e . (42) 

Thus, the action of B e on <j)p gives a function ipf +l \r,ee) oc Bi(f>^p(r), which satisfies 
the Schrodinger equation 

h e+1 ^=Ei^i £+1) - (43) 




012345 01234 



Figure 4: Left: The second excited eigenstate (p 2 ( r ) °f Ho (black curve) and the real 
(blue curve) and imaginary (red curve) parts of the corresponding Darboux-deformation 
■02 (r, e ) for e = 11 +«5. Right: The Born probability densities of 4>f\r) (black curve) 
and tp2 ] {r, e ). 



The straightforward calculation shows that this new function ips is square-integrable 
in Dy. Figure H] depicts the initial eigenfunction corresponding to i = 0, s = 2, as well 
as its Darboux-deformation ip^ (r, £0 = 11 +z5). The function <j>2 corresponds to the 
second excited state of H , accordingly it shows two nodes in Dy. However, although ^ 
corresponds to the second excited energy = 11 of hi, it is notable that this function 
is free of nodes. Thereby, the usual correspondence between the number of nodes and the 
level of excitation of a given energy eigenstate is missing in the new functions. Despite 
this fact, the normalized function can be put in correspondence with the Born's 

probability density. 

In the previous section we realize that the action of the canonical factorizing operators 
ae and a\ on a given physical state (j) s respectively reduces and increases in one unit the 
number of nodes (be, and b\, on the other hand, do not change the parameter s). In 
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contrast, nodes disappear in the complex Darboux-deformations of physical states, as 
established by the set of equations (|32|) . (1361) and (1421) . In this way, the parameter s in 
is nothing but the heritage of the number of nodes of the initial wave- function <f)*P. 
Another important consequence of the absence of this kind of zeros is that the new set 
of square- integrable functions } is not orthogonal. However, it is also possible to 

construct ladder operators which increase or decrease the energy of this system. 




Figure 5: Schematic representation of the intertwining operators Ae and Be, as they 
are applied on the energy levels of He and hi+\. The construction of the creation and 
annihilation operators of h,£+i is respectively depicted as the composites M\ +l = Be S e Ae 
and N e+1 = BeSeA £ . 



In analogy with (|27l) . let us analyze the intertwining relationships 

N t+1 (he+i - 4) = h l+l N l+l , M} +1 (h e +i + 4) = h e+1 M} +1 . (44) 
The above equations give place to the commutation rules 

[N £+1 , h t+l ] = AN i+1 , [M} +1 , h i+1 ] = -4M] +1 (45) 

which are the same as those fulfilled by H e , Sj and St in (EH]). Thus, M] +1 and N e+l 
should respectively play the role of the raising and lowering operators for the eigenvalues 
of he + \. From (jUJ), it is clear that Ni+xipg 1 ' is eigenfunction of he+i with eigenvalue 
E { s e) - 4 = Thus, N e+1 lowers the energy of the system in four units by annihilating 

one unit in the parameter s. Hence, we write V'i^i 1 '' °£ A^+l^s ■ m a similar way we 
get V's+i 1 '* Mj +l ipi i+1 \ Let us remember that S £ and Sj were constructed in terms of 
vertical displacements of the points in the (s x £)- plane. In the present case, we have 
to solve the algebraic equations (jUHlSD and Figure is of special utility. The simplest 
solution is given by the following products 

M\ +1 = B e Si Ae, N e+1 = B e S e A e , {M\ + ^ + N i+1 , (46) 

which can be verified by algebraic operations. Since Mj +1 and Ne+i are fourth order dif- 
ferential operators, they do not factorize the Hamiltonian he- The corresponding products 
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are fourth degree polynomials factorized as follows 



M} +1 N i+1 = (h t+1 -eg- A)(h t+1 + 2£- l)(h t+1 -21- 3)(h e+1 - eg) 

(47) 

N e+1 M\ +1 = (Vi -e e + 4)(h e+1 -2£+ l){h e+1 + 2£ + 3){h e+1 - eg) 

Finally, it is remarkable that, in contrast with Sj and Sg, the ladder operators Mj +1 and 
Ng +1 are not mutually adjoint (it is a heritage of the factorizing operators Ag and Bg). 
A similar situation gave place to the Distorted Heisenberg-Weyl Algebra reported in [24]. 
Research in this direction is in progress. 



4 Concluding Remarks 

Conventional factorization is useful to solve the eigenvalue problem of the spherical sym- 
metric harmonic oscillator. The corresponding Hilbert space 7i = L 2 (Dy) is spanned 
by the physical solutions which, in turn, are written in terms of orthonormal associated 
Laguerre polynomials. A new kind of factorizing operators has been introduced to express 
the radial oscillator Hamiltonian Hg as the product of two not mutually adjoint operators 
plus a complex constant. The reversed factorization produces a new Hamiltonian which is 
not Hermitian under the inner product of the associated Laguerre polynomials. However, 
the corresponding eigenvalue equation has been solved for the radial oscillator energy 
spectrum. Although the new solutions are not orthogonal, they are square-integrable in 
7i. Hence, the wave-functions so obtained can be put in correspondence with the Born 
probability density. 

The complex potentials of the non- Hermitian Hamiltonians reported here are such that 
their imaginary part is not a simple constant but a nontrivial function of the position. 
There is a well defined subset of Dy = [0, oo) in which this 'complex deformation' of 
the radial oscillator potential is relevant. A similar situation was reported for complex 
Coulomb-like [20, 21] and complex linear oscillator-like [19] interactions. Recent results 
show that, in a 'free particle background', this kind of potentials behaves as an optical 
device which both refracts and absorbs light waves [23] . Then, it is interesting to analyze 
the situation in a different background, namely in the radial oscillator one. On the 
other hand, in Section 3, it has been shown that the wave-functions of these complex 
potentials do not have a clear correspondence between the number of nodes and the level 
of energy excitation. However, our current studies show that this correspondence should 
be established by means of the real and imaginary parts of the new functions [33] (see 
also [34]). The improvements of these last results are going to be reported elsewhere. 
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